Real-space Manifestations of Bottlenecks in Turbulence Spectra 
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An energy-spectrum bottleneck, a bump in the turbulence spectrum between the inertial and 
dissipation ranges, is shown to occur in the non-turbulent, one-dimensional, hyperviscous Burg- 
ers equation and found to be the Fourier-space signature of oscillations in the real-space velocity, 
which are explained by boundary-layer-expansion techniques. Pseudospectral simulations are used 
to show that such oscillations occur in velocity correlation functions in one- and three-dimensional 
hyperviscous hydrodynamical equations that display genuine turbulence. 
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The energy spectrum E(k) characterizes the statistical 
distribution of kinetic energy among the wavenumbers k 
in homogeneous, isotropic, fluid turbulence in three di- 
mensions (3D). If k[ and kd denote, respectively, the 
wave-vector magnitudes related to the inverses of the 
lengths Li, at which energy is injected into the system, 
and rjd, where viscous dissipation becomes significant, 
then, in the inertial range k\ <C k <C kd, this spec- 
trum scales as E(k) oc k~ n ; the phenomenological theory 
(K41) of Kolmogorov which does not account for in- 
termittency 0], yields n = 5/3 for 3D fluid turbulence. 
In the far dissipation range k 3> kd, this spectrum falls off 
exponentially (upto algebraic prefactors) [3j. For values 
of k that lie in between inertial and far-dissipation ranges, 
a plot of the compensated energy spectrum (E(k) divided 
by its inertial-range form) versus k exhibits a gentle max- 
imum that is called a bottleneck 0, @| • Such bottlenecks 
have been seen in a variety of experiments @, @] and in di- 
rect numerical simulations (DNSs) of fluid turbulence @- 
[l(jl | . Phenomenological mechanisms have been suggested 
for the formation of bottlenecks (see, e.g., Refs. [!, [f|). 
A systematic theoretical study of the bottleneck phe- 
nomenon has been initiated in Ref. ll| by using the limit 
of high dissipativity a in hyperviscous hydrodynamical 
equations, which have a dissipation operator cx (— V 2 )". 
Hyperviscous dissipation, with moderate values of a, say 
between 2 and 4, is often used in DNSs in the hope of 
enhancing the inertial range of scales, but at the price of 
producing increasingly strong bottlenecks (see Ref. 11 1 
and references therein). 

In the first part of our study we develop a quantita- 
tive, analytical understanding of bottlenecks, for moder- 
ate values of cv, in the following hyperviscous generaliza- 
tion of the one-dimensional (ID) Burgers equation: 



d t u + ud x u = -v a k- 2a {-d 2 x ) a u + /fx, t); 



(1) 



here u{x,t) is the velocity at the point x and time t, 
v a > the hyperviscosity, k r a reference wavenumber, 
and / the driving force. It is well known that the ordi- 
nary (a = 1) Burgers equation, with / = 0, is integrable 



,12]; it is also easy to show that its energy spectrum has 
no bottleneck. By contrast, we show that, for any in- 
teger a > 1, the solution to the hyperviscous Burgers 
equation ([TJ, in the limit of small v a , displays an energy- 
spectrum bottleneck; for this it is crucial to examine the 
solution in real space, where we can use boundary-layer- 
type analysis, in the vicinities of shocks, to uncover oscil- 
lations in the velocity profile. We obtain this result both 
for the unforced, hyperviscous Burgers equation and for 
its variant (DHB) with deterministic, time-independent, 
large-scale forcing. We validate our DHB solutions with a 
pseudospectral DNS. Note that these solutions are time- 
independent and not turbulent; however, the key qualita- 
tive feature of real-space oscillations in the velocity pro- 
file does carry over to oscillations in velocity correlation 
functions in one- and three-dimensional hyperviscous hy- 
drodynamical equations that display genuine turbulence. 
We show this in the second part of our study by using 
DNS. This association of bottlenecks and oscillations in 
velocity correlation functions has not been made so far. 
It is akin to the association of peaks in the static struc- 
ture factor S(k), of a liquid in equilibrium, with damped 
oscillations in the radial distribution function g(r) [13j. 



The simplest model presented here, which displays a 
bottleneck amenable to analytical study, is the DHB 
equation (p} with integer a > 1, time- independent force 
/ = sin £, and u(x,t — 0) = 0. The velocity eventu- 
ally goes to a steady state, which is a solution of the 
ordinary differential equation (ODE) that is obtained by 
dropping the time-derivative term in ([lj. When a ^ 1, 
this nonlinear ODE is not integrable, but its limit as 
v a — > is the same as for ordinary dissipation, namely, 
it has a shock at x = 7r, where the solution jumps from 
tt_ = +2 to u + = —2. For small but finite v a , the 
shock is broadened and its structure can be analyzed 
by a boundary-layer technique using the stretched spa- 
tial variable X = (x — ix)jv^ , with (3 — 2 a- \ > an< ^ ex ~ 
panding the boundary-layer velocity u BL in powers of v a : 
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m bl (X) = £jl z4u,(X). To leading order (j = 0) 

M*)-<- iyM £*» «•<*»» -=*■ < 2 > 

For a = 1, we obtain the standard profile uq = 
— 2tanhJT. For a > 1, Eq. ([2j cannot be solved an- 
alytically. However, for large X the equation can be 
linearized because uo is close to its asymptotic con- 
stant value. For example, for large negative X, we set 
Mo = 2 + w, discard the quadratic term in w, and obtain, 
after integrating once, (— l) a+1 d 2a ~ 1 w / dX 2a ~ x = 2w. 
This constant-coefficient ODE has solutions of the form 
fiexp(n a X), where fi is arbitrary and the "eigenvalue" 
n a is any of the (2a - l) th roots of (-1) Q+1 2, i.e., 
for even a, n a = 2 /3 exp(i(2n + 1)13%) and for odd a, 
K a = 2 13 exp(z2n/?7r), with n = 0, 1, (2a - 2). Only the 
eigenvalues that have a positive real part are acceptable, 
because w should vanish at — oo. If all the modes with 
such eigenvalues are actually present then, for X — > — oo, 
the solution to ((2]) tends to +2 in an oscillatory fashion 
and it is dominated by the mode n* (and its complex 
conjugate), which has the smallest positive real part. In 
terms of the unstretched coordinates, this means that, in 
the neighborhood of the shock, the solution for even a 
displays damped oscillations with wavelength 

A* = 2™f [2" sin[(2n* + 1)/3tt)]] -1 (3) 
and with an e-folding rate 

tf* h = 2^cos[(2n* + l)/?7r)]. (4) 

The case of odd a is handled mutatis mutandis. 

Such damped oscillations imply the presence of a pair 
of complex k poles in wave-number space, whose signa- 
ture, for real k, is a Lorentzian. This can be a bump or a 
trough, near wave number 27r/A Q , with width ~ K a and 
amplitude ~ K~ . We present below a semi- numerical 
analysis to show that the solution of the DHB yields a 
bottleneck (bump). 

The theory presented here can be applied to a wide 
class of problems ranging from the unforced, hypervis- 
cous Burgers equation to the case of arbitrary, large- 
scale, time-dependent or time-independent forcing as 
long as it does not modify the internal structure of 
shocks. Note also that such a linearized theory gives 
no prediction for the a-dependence of the amplitude of 
the bottleneck. Furthermore, when more than one mode 
(or pair of complex-conjugate modes) with a positive real 
part is present (i.e., for a > 3), linear theory does not 
tell us if the mode(s) with the smallest positive real part 
is (are) actually excited. Such issues require a global 
analysis of the boundary-layer equation ([2j and not just 
of its large-X, linearized version. Except for the very 
standard case of ordinary dissipation, we do not know 
much about the properties of the solution to Eq. ©. 



If we demand that Uq(X — 0) = (which can always 
be achieved by a suitable translation), is the solution 
unique? By using a numerical, shooting method, we 
obtain evidence that, for a = 2, there is a unique so- 
lution that has u' (X)\ x =o = -2.121530817618... and 
u'q(X)\x=o = 0. We can also obtain the value of this first 
derivative at the origin with ~ 10% accuracy by assum- 
ing that the solution has singularities on the imaginary 
axis at X = Z± = ±i A (a Painleve-type argument indi- 
cates that, near such a singularity, to the leading order, 
uq(X) w 120 /(X — Z+)). The vanishing of the second 
derivative implies that this unique solution is odd in the 
X variable. 

Direct numerical integration of the boundary-layer 
equation ([2} is a greater challenge than the full DHB 
equation because, for the latter, we can take advantage 
of periodic boundary conditions. By using the value of 
Uo(0), obtained by the shooting method discussed above, 
and u'q(0) = 0, we solve the third-order, boundary-layer 
equation for a = 2 numerically. We find Xc — 1.15, the 
value of X at which uo(X) first crosses the —2 asymptote. 
Next, we calculate u' (Xc) by using the Taylor expan- 
sion u a (X) ~ u' (0)X + u "(0)X 3 /3! + w ""(0)X 5 /5! + ... 
along with the known values of u' (from the shoot- 
ing method), and <'(0) = 2 and <"'(0) = -(u (0)) 2 
(from Eq. @). The linear theory suggests Uo(X) = 
-2 + Ae- K2( > x - Xc '> sin2ir(X - X c )/\2 for X > X c ; 
thence we obtain u' a (Xc) in terms of A and A2. By using 
the values of u' (Xc) (from the Taylor expansion above) 
and A2 (cf. Eq. ([3])), we obtain A ~ —0.983, which is 
within 1.7% of the value of A (~ —0.966) that we get 
from the solution of the boundary-layer equation. 

We now address the question of whether the Fourier- 
space manifestation of these oscillations is a bump or a 
trough. The Fourier transform j(k) of the real and even 
function —u'(X) is real and even [l9j]; and j(k) is the 
square root of the compensated energy spectrum. The 
rising of the compensated energy spectrum, in the inter- 
mediate regime, between the fiat region near k = and 
the exponential decay at large k, is equivalent to j"(k) 
being positive; and 7"(fc)| fe ^ = l/2n J °° dXX 2 u'(X). 
To solve for Y'(k), we use u'(X) either from a numerical 
solution of Eq. ([5]) or from the linear theory above; we 
then perform a numerical integration over X\ we obtain 
good agreement (~ 9%) between the results of both these 
methods; and, indeed, we find that "f"(k) is positive, so 
the spectrum has a bottleneck. 

We turn now to a pseudospectral DNS of the DHB 
equation ([1} with a = 2,4,8, and 16, a 2/3 rule for 
dealiasing, and a fourth-order, Runge-Kutta method for 
time-integration. The reference wavenumber k r = 100; 
the number of collocation points N = 2 14 , the time 
step St = 10~ 4 , and the hyperviscosity coefficients are 
v 2 = 5 x 10~ 3 , v 4 = 5 x 10~ 8 , v% = 5 x 10~ 14 , and 
vie = Kr 20 . 

The steady-state, compensated energy spectra E% = 
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FIG. 1. (Color online) Plots for the DHB Eq. fl}: (a) Log-log plots of the compensated energy spectrum k 2 E(k) versus k 
for a = 2 (black squares), a = 4 (blue filled-circles) , a — 8 (cyan diamonds), and a — 16 (red hexagons), (b) Plots of the 
steady-state solution u(x) for the same values of a as in (a); these are indistinguishable from u°(x) (thick magenta line) away 
from x — 7r; inset: plots of u d (x) versus x around x = ir. (c) Semilog plot of Aw (red filled circles) versus ir — x; the black line 
is the fit. 



k 2 E(k) (Fig. (la)) show clear bottlenecks; the height 
of the bottleneck peak increases with a; but k^ , the 
wavenumber of this peak, decreases as we increase a. 
We now investigate the real-space manifestation of this 
bottleneck. In Fig. (lb) we plot the steady-state so- 
lution of the DHB equation; this shows that our nu- 
merical solution agrees with the outer solution u°{x) = 
2sgn(7r — x)sm(x/2) away from the shock. However, 
in a thin boundary layer around the shock at x = ir, 
there are conspicuous oscillations that become promi- 
nent when we plot the difference between the solution 
of the DHB equation and the outer solution, u d (x) = 
u(x) —u°(x), versus x [Fig. (lb), inset]. The character- 
istic wavelength of these oscillations is X a ; for the rep- 
resentative case a = 16, we find Ai6 = 0.0122. Simi- 
larly, from plots of the compensated spectra [Fig. (la)] 
we obtain Ai6 = 0.0121. Furthermore, the theoretical 
prediction for the wavelength of these oscillations (cf. 
Eq.©) yields Ajg = 0.0120. The bottleneck has a fi- 
nite width because of the decaying envelope of the oscil- 
lations in real space as we move away from the shock (cf. 
Eq.©). We obtain the amplitude A a of these oscilla- 
tions [Fig. (lb)] as a function of (x — ir) and find, numer- 
ically, that A a ~ exp[K a (x — ir)] as shown in Fig. (lc). 
For a = 16 we obtain, from our DNS, an e-folding rate 
Ki 6 ~ 26.61, whereas our theoretical prediction in Eq.© 
yields K{g ~ 26.54. Thus, we find excellent agreement 
between our theoretical predictions (Eqs.© and ©) and 
our numerical results for both the wavelength of the os- 
cillations and the e-folding rate. 

Our results carry over to the stochastically forced hy- 
perviscous Burgers equation (SHB) and to the 3D hyper- 
viscous Navier-Stokes (HNS) equations, both of which 
exhibit turbulence. Let us first examine bottlenecks in 
the SHB equation jl4j |. with a white-in-time, Gaussian 
random force with zero mean, an ultraviolet cutoff at 
N/8, and a spectrum ~ k . The velocity field for the 
SHB shows shocks at various length scales; and the re- 
sulting energy spectrum shows an inertial-range scaling 
E{k) ~ fc -5 / 3 . In Fig. (2a) we give a representative plot 



of the compensated energy spectrum k 5 ^ 3 E(k), for a = 8; 
this shows such an inertial range followed by a promi- 
nent bottleneck that peaks at a wavenumber kf ~ 890. 
We measure the correlation function (u(x)u(x + I)) , aver- 
aged over statistically independent configurations, which 
show oscillations [Fig. (2a), inset]; these are the real- 
space manifestations of this bottleneck. The wavelength 
of these oscillations is ~ 0.00706; and the correspond- 
ing wavenumber is ~ 889.97, in agreement with the 
wavenumber at which the bottleneck shows a peak in 
Fig. (2a). 

The 3D hyperviscous Navier-Stokes (HNS) equation 
for an incompressible velocity field u(x, t) is 

+ u • Vu(x, t) = - Vp - v a (- V 2 ) ° u(x, t) + f (x, t); 

V • u = 0. (5) 

We integrate it by a pseudospectral method with a 2/3 
dealiasing rule, an Adams-Bashforth scheme for time 
marching, 512 3 collocation points, a = 4, and — 
10~ 14 . We force the 3D HNS equation to a statisti- 
cally steady state by using the constant-energy-injection 
method described in Ref. In Fig. (2b), we show a 

representative plot of the compensated energy spectrum 
E c (k) = k 5 ^ 3 E(k); this shows a bottleneck between the 
inertial and dissipation ranges. The correlation function 
D(l) = (u(x) • u(x + 1)), averaged over five configura- 
tions that are separated from each other by 677, where 
17 is the integral-scale eddy turnover time, shows gentle 
oscillations [Fig. (2c)], which are the real-space manifes- 
tations of this bottleneck. These oscillations can be seen 
clearly in D°{l) [Fig. (2c), inset], which is obtained by 
subtracting the linear, decaying trend from D(l). The 
wavelength of these oscillations is ~ 0.1665 and the cor- 
responding wavenumber is ~ 37.7, in agreement with the 
wavenumber at which the bottleneck shows a peak in 
Fig. (2b). 

We have provided a theoretical explanation for energy- 
spectra bottlenecks in the DHB equation by combin- 
ing analytical and numerical studies. These bottle- 



4 




FIG. 2. (Color online) (a) Log-log plot of the compensated energy spectrum k 5/l3 E(k) versus k for the stochastically forced 
hyperviscous Burgers (SHB) equation with a — 8 and (inset) a plot of its correlation function (u(x)u(x + l)) showing oscillations 
of wavelength A| HB , which is inversely related to the wavenumber of the bottleneck in the energy spectrum, (b) The compensated 
energy spectrum k 5 ' 3 E(k) for the 3D HNS equation (a — 4) with a bottleneck peak at wavenumber K™ s = 40. (c) A plot 
of the correlation function D(l) versus I for the 3D HNS equation; inset: oscillations in a plot of the function D°(l), which we 
obtain by subtracting the linear, decaying trend from D(l). 



necks appear as a natural consequence of oscillations 
in the velocity profiles in the vicinity of a shock. Ear- 
lier studies lg, [l7| have seen such oscillations in the 
DHB case but have not associated them with bottle- 
necks in energy spectra. Furthermore, we have shown 
that energy-spectra bottlenecks in the SHB and the 3D 
HNS equations, which exhibit turbulence, are associated 
with damped oscillations in real-space velocity correla- 
tion functions. This association has not been made hith- 
erto, even though there have been attempts to explore 
real-space manifestations of bottlenecks in energy spectra 
10]. Our work confirms that the larger the dissipati vity 
a, the more pronounced is the bottleneck 0, [HI Il8| . 
Thus, by increasing a, we have uncovered clearly the 
real-space manifestations of these bottlenecks. Energy 
spectra for homogeneous isotropic turbulence in the 3D 
NS equation (a = 1) show a mild bottleneck @, [HI; we 
expect, therefore, that there should be weak oscillations 
in real-space velocity correlation functions. The detec- 
tion of such weak oscillations is an important challenge 
for experiments and DNS. 
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